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Abstract 

This paper is devoted to general balance laws (with a possibly non local 
source term) with a non-characteristic boundary. Basic well posedness 
results are obtained, trying to provide sharp estimates. In particular, 
bounds tend to blow up as the boundary tends to be characteristic. 
New uniqueness results for the solutions to conservation and/or balance 
laws with boundary are also provided. 
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1 Introduction 

This paper is concerned with initial boundary value problems (IBVP) for 
systems of balance laws of the form 

dtu + d,f{u) = G{u) X > -fit) 
6(n(t,7(t))) =5(0 t > (1.1) 
u{0,x)=Uo{x) X > 7(0) 

where / is smooth, Df is strictly hyperbolic, Uq is the initial datum and 
G is a possibly non-local source term. The boundary 7 is assumed non 
characteristic, i.e. £ characteristics point outwards and n — i inwards. The 
role of b is that of letting n — £ component of u be assigned by the boundary 
data g. Above and in what follows, we assume that all BV functions are 
right continuous. 

Systems belonging to this class were already considered in the literature. 
See, for instance, [3 18] for the case with a non local source but no boundary 
and [12] for the case of a Temple type /. 
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Examples of physical models that fit into this class are found, besides 
in the cited references, also in [11]. There, a model describing the flow of a 
fluid in a simple pipeline is based on a system essentially of the form (II. ip . 

As is well known, preliminary to the study of (jl.ip . is that of the purely 
convective system 

' dtu + d^f{u) = X > 7(t) 

< 6(i.(t,7(t))) =g{t) t > (1.2) 
u{0,x)=Uo{x) X > 7(0) 

considered, for instance, in [U [21 [121 13 [S]. Below, we provide results 
on ()1.2p that are not contained in these papers. In particular, the present 
estimates explicitly blow up as the boundary tends to be characteristic. The 
choice of the Glimm type functionals on which most of the proof relies is 
here simplified, compare for instance (j4.1ip below with [15\ (2.10)-(2.13)] 
and gUD-gTSI) with [151 (3.5)-(3.10)]. 

In the homogeneous case (|1.2p , we also provide a uniqueness result that 
has no analogue in the case of Cauchy problems with no boundary. Indeed, 
let u solve (II. 2p and assume a second boundary 7 is given, such that 7 > 7. 
Along 7 assign the trace of u as boundary data, i.e. let g{t) = b (^u {t, ^{t)^^. 
Then, the solution to 

' dtu + d^f{u) = X > 7(t) 

< 6(n(t,7(t))) =m t > (1.3) 
u{0,x) = Uo{x) X > 7(0) 

coincides with the restriction of u to x > 7(t), see Proposition 12.41 We show 
that an analogous result may not hold in the case of (jl.ip . see (j3.ip . 

Besides, we also provide a Lipschitz estimate on the process generated 
by (|1.2p that contains also a second order part on a generic perturbation, 
see[^ in Theorem l2.2l This technical estimate, already known in less general 
situations, played a key role in several other results, see for instance [71 
Proposition 3.10] and [21 Remark 4.1]. 

All what we obtain in the case of ()1.2p is used in the proof of the results 
on p.ip . In particular, for both systems, we provide bounds on the total 
variation of time like curves. These estimates are optimal in the sense that 
they blow up as the boundary tends to be characteristic, see propositions 12.31 
and 13.31 

The next section is devoted to the homogeneous problem (jl.2p . while 
Section [3] presents the results related to (II. ip . The proofs are deferred to 
the last two sections. 
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2 The Purely Convective IB VP 



On system (|1.2p we require the following conditions: 

(/) /: 17 — > M" is smooth, with Q C R" being open, such that Df{u) 
is strictly hyperbolic for all u € fi, each characteristic field is either 
genuinely nonlinear or linearly degenerate. 

Without loss of generality, we may assume that € and for all u in Q, 
Df(u) admits n real distinct eigenvalues Ai(n), . . . , A„(u), ordered so that 
\i-i{u) < Xi{u) for all n, with right eigenvectors ri(n), . . . ,rn{u). 

(7) 7 G C°'^(]R+;R) and, for a fixed positive c, Xi{u)+c < ^{t) < Xe+i{u)- 
c for a fixed £ G {1, . . . , n — 1} and for all ti G il. 

(b) b G C^{n;M."''^) is such that 6(0) = and 

det [DbiO) n+i (0) DbiO) ri+2 (0) • • • D5(0) r„ (0)] / . 

For notational simplicity, we say below that a curve 7 is £-non- characteristic 
if 7 G C°'-'-(IR+; M), and for a fixed positive c, for all u G fi, A£(u) + c < 
i{t) < A^+i(ti) — c. This notion is more restrictive than that of a non- 
resonant curve, see [131 Chapter 14]. 
We define below the domain 



and extend to [0, +oo[ x IR any function defined on to vanish outside D-y. 

We slightly modify the definition given in [IT] of solution to ()1.2p in the 
non characteristic case, see also [H [21 [T^] and [HI Definition 2.1]. Indeed, 
here we require the boundary condition to be satisfied by the solution only 
almost everywhere. This softening allows for a simpler proof without any 
substantial change, since we provide below a full characterization of this 
solution, seed]), [2]) with a; = and [3]) in Theorem 12. 2[ 

Definition 2.1 Let T > 0. ^ map u = u{t,x) is a solution to if 

1. u e CO ([0,r];Li(M;M")) with u{t,x) G n for a.e. {t,x) G and 
u{t, x) = otherwise; 

2. n(0, x) = Uo{x) for a.e. x > 7(0) and lim b(u{t,x)) = g{t) for 

x— >o+ 

a.e. t> 0; 

3. for X > 7(t), u is a weak entropy solution to dfU + dxf{u) = 0. 



D, 
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Theorem 2.2 Let the system satisfy (f), (h), Assume also that 
g G BV(M+; M"'^^) has sufficiently small total variation. Then, there exists 
a family of closed domains 

Pt C |u G (L^ n BV) (M; Q) : u{x) = for all x < 7(t)| 

defined for all t > and containing all functions with sufficiently small 
total variation that vanish to the left ofjlt), a constant L > and a process 

P{t,to):Vt^ ^Vt^+t, for allto,t>0, 

such that 

1) for all to > and u G Vf^, P{0,to)u = u while for all t,s,to > and 
ueVt^, P{t + s,to)u = P {t, to + s)o P{s, to)u; 

2) let uj he an L-*- function with small total variation, if (Pj'Dt) are the 
process and the domain corresponding to the boundary j{t) and boundary 
data g{t), then, for any u G Pj^, v G 'T't'^, we have the following Lipschitz 
estimate with a second order error term accounting for to: 

\\P{t,to)u-P{t',t'o)v-i0\\-^^^ 

< L ■ i^\\u - V - uj\\^i + \t -t'\ + \to -t'^\ 

j-to+t 

+ \\9{t) - g{T)\\dT + sup |7('^)-7(^)| 

Jto Te[to,t] 

+t-TVH I; 



3) the tangent vector to P in the sense of 15^ Section 5] is the map F defined 
at l^.gp , i.e. for all to > and u G T>t^ 

\\m.J^\\F{t,to)u - P{t,to)u\\^^ =0; 

4) for all Uo G Vq, the map u{t,x) = (P(t,0)iio) (x) defined for t G [0, T] 
and {t,x) G B^, solves U.2\) in the sense of Definition \2. 1\ 

P is uniquely characterized byUl), with w = and\^. 

The conditions [7p-[3) constitute what is the natural generalization to the 
present case of the definition of Standard Riemann Semigroup., see Defi- 
nition 9.1]. 

Remark that, in general, the Lipschitz constant L blows up as the the 
boundary tends to become characteristic, i.e. as c — > 0, see (|4.18p and the 
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next proposition. Indeed, in the proof of Theorem 12.21 we prove also the 
following result on the regularity of the solutions to ()1.2p along non charac- 
teristic curves. 



Proposition 2.3 Fix a positive T. Let the system U.^) satisfy the assump- 
tions of Theorem \2.2\ and callu the solution to constructed therein. Let 
ro,ri be i -non- characteristic curves, for I £ {!,..., n — 1}. Then, there 
exists a constant IC > independent from T,Uo,g such that 

u(t,ro(t)) -n(t,ri(t)) (it < -(TVK) + TV(5)) ||ri-ro||co([o,T])- 



A uniqueness property proved in Section S] is the following. 

Proposition 2.4 Let the system 
satisfy the same assumptions of Theo- j. 7(^) 
rem \2.S\ and call u the solution to 
constructed therein. Let 7 E C'''-'-(M+; M) 
he any £ -non- characteristic curve satis- 
fying 7(t) > 7(t) for all t > 0. De- 
fine g{t) = b{u{t,^{t)+). Then, / TO) 
also satisfies the assumptions on Theo- 
rem \2.2\ and the solution u constructed by 
this Theorem satisfies 



lit) 



u{t, x) 



u{t, x) 
n(t, x) 



u{t, x) = u{t, x) 
for all X > 7(t) and t >0. 

3 The IBVP with General Source Term 

To deal with the source term, for all positive 6, define 

= j-u G L^(R; Q) : TV(n) < 5^ . 

We add the following assumption on the source term of system 

L-'-(IR,M") is such that for suitable positive 



(G) For a positive 60, G: Us^ 
Li,L2 



\G{u) — G{w)\\-^i < Li ■ \\u — w\ 
TY{Giu)) < L2. 



The natural extension of Definition 12.11 to the present case is the following. 
Definition 3.1 Let T > 0. A map u = u{t,x) is a solution to if 
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-Z. n G C° ([0,T];Li(M;M")) with u{t,x) € n for a.e. {t,x) £ 
u{t, x) = otherwise; 



and 



2. u{0,x) = Uo{x) for a.e. x > 7(0) and lim b(u{t,x)) = g{t) for 

3:^0+ 

a.e. t> 0; 

3. for X > 7(t), u is a weak entropy solution to dfU + dxf{u) = G{u). 

With this notation, we may now state the extension of Theorem 12.21 to the 
present non homogeneous case. 

Theorem 3.2 Let system satisfy (f), (G), (b), ('y). Assume also 

that g € BV(M'^;R"~^) has sufficiently small total variation. Then, there 
exist positive 5, L, T, domains T>t, for t € [0, T] and maps 

P{t,to): Vt„^Vt„+t 
for to,to + t € [0, T], such that 
i) {ueUs: u{x) =Oforx< -fit)}; 

a) for all to, ti,t2 with to G [0, T[, ti € [0, T — to[ and t2 E [0,T — to — ti], 
P{t2,to + ti) o P{ti,to) = P{ti + t2, to) and P(0, to) = Id; 

Hi) if{P,'Dt) are the process and the domains corresponding to the boundary 
7(t) and boundary data g{t), satisfying the same assumptions above, 
then, forto,t'o G [0,T[, t G [0,T-to] andt' £ [0,T-t'o], for all u G 
u G Vt' 



P{t,to)u - P{t' ,t'o)u 



< L ■ {\\u — u\ 



1 + ll-"! 



\t, t, \ \ tQ t/Q 



to+t I 

+ / \\9{t) - g{r)\\dT + sup |7(t) - 7(t)| >; 

'to T&[to,t] I 



iv) for all to G [0,T[, t G [0,T — to], u G Vt^ define 



F{t, to)u = P{t, to)u + t G{u) X 



[■y{to+t),+<x>[ 



then 



lim — 

t^o+ t 



P{t,to)u- F{t,to)u 



0: 



v) for all Uo G Vq, the map u{t,x) = \^P{t,0)uoj (x) defined for t G [0, T] 
and {t,x) G B^, solves U.l\) in the sense of Definition \3. 1\ 
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The process P is uniquely characterized bu\i^). \m\) and \vt^) . 



Again, as remarked after Theorem 12.21 the Lipschitz constant in general 
blows up as c — > 0. The proof of this result is deferred to Section [SJ it 
heavily relies on Theorem 12.21 Remark that it is possible to extend to the 
non homogeneous case also Proposition [ 



Proposition 3.3 Let system fHiP satisfy the same assumptions of Theo- 
rem \3.S\ and call u the solution to constructed therein. Let ro,ri be 

i -non- characteristic curves, for I G {1, . . . , n — 1}. Then, for all Uq and g, 
there exists a constant /C > such that 



u{t,ro{t)) -u{t,ri{t)) 



dt< — \\Ti 
c 



■ 0|lco([0,T])- 



Remark 3.4 Proposition \3.3\ implies also that, ifV is any I -non- characteristic 
curve, then the map x (^P^t, 0)^^ (x) is continuous in x = T{t) for almost 

all t G [0,T]. Indeed, denote u{t,x) = (P{t,0)u] (x) and compute 



T 







u{t,T{t)-) -u{t,T{t)) 



dt 



1 



u{t,r{t)-x) -u{t,T(t)) 



dt dx 



= lim — 

£ JQ JQ 

/c , 1 r , 

< — lim - / xdx = 

C e^O £ Jo 

Contrary to Proposition 12.31 the uniqueness result of Proposition [27 
may not be extended to the present non homogeneous case, due to the non 
local nature of the source term here considered. Indeed, let 
J 7(0 ^^*K{t,x) 

dtu + d^u= Jq u{t, ^) 
u{t,0) = 




(3.1) 



u{0, x) 



It is immediate to verify that the assumptions of Theorem 13.21 hold. The 
solution u, shown above, is non zero in the delimited area above and, in 
particular, for t £ [0, 1] and x £ [3, 4] but it vanishes for t £ [0, 1] and x = 1. 
Therefore, with the same notation of Proposition 12.31 letting 'y{t) = and 
7(t) = 2 we have g{t) = for t G [0, 1]. Problem (II. 3p thus admits, in the 
present case, only the trivial solution u = 0, contradicting what would be 
the analog of Proposition 12.31 in the non homogeneous case. 
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4 Proofs Related to Section [2] 

Below, C denotes a positive constant dependent only on /, G and b whose 
precise value is not relevant. 

This section is devoted to the homogeneous initial boundary value prob- 
lem (ll.2p and proves Theorem 12.21 Our general reference on the theory of 
conservation laws is [6]. 

Let a —5- Rj{a){u), respectively a — > Sj{a){u), be the j-rarefaction curve, 
respectively the j-shock curve, exiting u. If the j-th field is linearly degen- 
erate, then the parameter a above is the arc-length. In the genuinely non- 
linear case, see [U, Definition 5.2], we choose a so that (see [U formula (5,37) 
and Remark 5.4]) 



{R,{a){u)) = 1 and ^ (S,(a)(n)) = 1 . (4.1) 



Introduce the j-Lax curve 



Rj{a){u) if cr > 
Sj{a){u) if < 



and for a = (cji, . . . , (j„), define the map 

*(cr) = tpnicFn) O ...O ■0l((Ti) . 

By (/), see [U Paragraph 5.3], given any two states G 17 sufficiently 

close to 0, there exists a map E such that 

a = E{u^ ,u^) if and only if u+ = *(<t)(u") . (4.2) 

Similarly, let the map S and the vector q = (gi, . . . , g„) be defined by 

u+ = S{q){u-) and 5(g) = 5„(g„) o . . . o 5i(gi) , (4.3) 

i.e. S is the gluing of the Rankine - Hugoniot curves. 

We first consider the non characteristic Riemann problem at the bound- 
ary 

^ dtu + d^J{u) = Q X > 7(t) 

b(u{t,-i{t))^=go t > (4.4) 
u{0, x) = Uo X > 

where go G M"~^ and G are constants and 7 satisfies (7). Then, a 
solution to (14. 4p is constructed as in , see also [U El |12] . 
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Lemma 4.1 Let (f), (^) and (h) hold. Ifuo,go o-Te sufficiently small, then 



there exists unique E7,E? of class and states u'^,u'^ such that 



{cTg+i, . . . ,(T„) = E^{uo,go) 

{qe+i,---,qn) = El{uo,go) 



b{u'^) = Qo and 

IpnicTn) o . . . o V^4.i(cr^+i)(u<^) = Uo , 

b{u'^) = go and 

Sn{qn) O ...O S'^+i(g£+i)(M«) = Uo ■ 



Proof. We prove this statement only for the Lax curves, the results for 
the shock curves is proved similarly. Let cxj — > ■i}ji{ai){u) be the inverse Lax 
curve, i.e. 

a^^.(a)(u) - I "^J-^^^^^) if a > 
a^y;j[a)[u) - j a < 0. 

The choice (14. ID of the parameters implies that ipi{—cri) o %l)i[(ji){u) = u for 
all small u and ctj. Define the function 

G(cr^+i, . . .,an,go,Uo) = b {ipi+ii-ai+i) o ... o ^„(-cr„)(Uo)) - go- 
By (b), G satisfies G(0,0,0) = and 

det ^^^)G(0, 0, 0) = (-1)-^ det [Db{Q) r,+i (0) • • • 1)6(0) r„(0)] / . 

The Implicit Function Theorem guarantees the existence of a map E^ = 
E^{uo,go) with the required properties, if (cr^+i, . . . , cxn) = E^{uo,go) and 

= tpe+li-CFi+i) o . . . o 1pni-Crn)iUo). □ 

The notation introduced above allows the definition of a local fiow tan- 
gent to the process generated by ()1.2p . Fix to > and u £ T>t^, define 
go = g{to+), Uo = u (7(^0)+) and u'^ as in Lemma HTTl Let 

I U[X) it X > 7(^0) 

Call S the Standard Riemann Semigroup generated by /, see [HI Defini- 
tion 9.1]. Finally, for t > 0, define the tangent vector, see [5l Section 5], 

We record here the following interaction estimates, see Figure [H 
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Lemma 4.2 Let (f), (^) and (h) hold. If the following relations hold 
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b{u-) 

b{u+) 



Ur 



IpnicTn) o ■ ■ ■ oipli(Tl)iu ) 



then, we have the estimate 



- g' 



i=l 



Analogously, for the shock curves, if uj is a small vector satisfying 
g = Ku), 



g = h{v), 
then, we have the estimate 



and V + ijO = Sn{qn) o . . . o Si{qi){u) 



i=e+i 



i=l 




ig,g) 



{u,v) 



Figure 1: Interactions at the boundary 

Proof. By Lemma (ct^+i, . . . , = E^{ur,g'^) and (a^+i, . . . , cr„) = 
(^Ur, b [Tpe{(7i) o • • • o 'ipi{ai){u~))^ . Therefore, the Lipschitz continuity 
of imphes: 



i=e+i 



XI \di-ai\ < ci b{u^) -b(^Tpe{ae)o---oilji{ai){u~fj 



+ 



9 -9 



9 -g 



1=1 



Concerning the shock curves, by Lemma 14.11 we can write 

{qi+i,...,qn) = E^(^v + u;,b{Si{qi)o---oSi{qi){u))^ . 
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Since {v, g) = and 6(u) = g, the Lipschitz continuity of and b implies 

n 

< C El(v + u;,b{Se{qe)o---oSi{qi){u)))-El{v,g) 



< C (^\\oj\\ + \\b (Seiqe) o . . . o Si{qi){u)) - b{u) 



+ \\9 - g\ 



< 



c ^\qi\ + Ikll + 



\9 - 9\ 



completing the proof. 



□ 



Remark 4.3 We need below the first statement of Lemma \4.2\ in the partic- 
ular case of only one incident wave, i.e. Ur = ipi{(Ti){u'^) for some 1 < i < i. 

We follow the nowadays classical wave front tracking algorithm, see [U 
El El [12] , to construct solutions to the homogeneous boundary value prob- 
lem ([12]). Let M G (]7(i),+oo[,M") be piecewise constant with finitely 
many jumps and assume that TV(n) is sufficiently small. Call J{u) the 
finite set of points where u has a jump. Let crx,i be the strength of the i-th 
wave in the solution of the Riemann problem for 



dtu + dxf{u) = 



(4.7) 



with data u{x—) and u{x+), i.e. {ctx^i, . . . , cyx.n) = E (n(x— ), u(x+)) . Obvi- 
ously, if X ^ J(n) then ax,i = 0, for alH = 1, . . . , n. In x = 7(t) define 



CJ, 



7(i),l' 



,0". 



7(t),ra 



,0,i?,-(n(7(t)+),<7(t))) • (4.8) 



Then, consider the Glimm functionals and potentials 



a:>7(i) «=1 x>^{t) i=l+l 

Ttiu) = Vtiu) + H2Qt{u) + HiTY{g,[t,+<^[} 



(4.9) 



the set A of approaching waves being defined as usual, see [6] and the con- 
stant K,Hi,H2 to be defined later. As in [S], using Lemma [4.21 the Glimm 
functional Y can be extended in a lower semicontinuous way to all functions 
with small total variation in L-'-(IR;]R"') that vanish for x < "y{t). On the 
contrary, the interaction potential Q alone does not admit a lower semicon- 
tinuous extension, due to the presence of the boundary. 
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We now construct e-approximate solutions to (II. 2p by means of the classi- 
cal wave front tracking technique, see [6] or [Tl|2] for the case with boundary. 

Let e > be fixed and approximate the initial and boundary data in ()1.2p 
by means of piecewise constant functions and such that (see [9, for- 
mula (3.1)]) 



^ ^ (4.10) 
|To«)-To(no)| <e, TV(/,]t,+oo[) <TV(5,]t,+oo[) for t > . 

To proceed beyond time i = 0, we construct an approximate solution to ()1.2p 
by means of the Accurate and Simplified Riemann solvers, see [6l Para- 
graph 7.2]. Introduce the threshold parameter p > to distinguish which 
Riemann solver is used at any interaction in x > ^(t). Whenever an interac- 
tion occurs at (t,x) with x > 7(t), proceed exactly as in [6[ Paragraph 7.2]. 
Recall that the former solver splits new rarefaction waves in fans of wavelets 
having size at most e, while the latter yields nonphysical waves. These waves 
are assigned to a fictitious n+l-th family and their strength is the Euclidean 
distance between the states on their sides. 

At any interaction involving the boundary, i.e. when a wave hits the 
boundary as well as when the approximated boundary data changes (see Fig- 
ure [H left and center), we use the Accurate solver, independently from the 
size of the interaction. As usual, rarefaction waves are not further split at 
interactions. 

Along an e-approximate solution, for suitable constants K,Hi,H2 all 
greater than 1, introduce the linear and quadratic potentials and the Glimm 
functional: 



t n+l 



x>'y{t) i=l x>y{t) 1=1+1 



Vg%t) = TV(/;[t,+(X)[) (4-^^) 

T%t) = V'{t)+H,Vg%t) + H2Q%t). 

The potentials just defined differ from the ones defined in ()4.9|) because the 
non-physical waves are accounted for in a different way. They coincide at 
t = because of the absence, at that time, of non-physical waves. They 
differ of a quantity proportional to the total size of non physical waves for 
t > 0. 

As usual, changing a little the velocities of the waves, we may assume 
that no more than two waves a' ,a" collide at any interaction point it,x). 
When X > 7(t), the usual interaction estimates yield, for a constant C > 
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dependent only on / and 17, 



AV%t) < CK\a'a"\ AF/(t) = 



as soon as CK < -ff2/4 and 5o is sufficiently small. 

When a wave a hits the boundary, Lemma 14.21 implies that 



AV%t) < {C-K)\a\ AVg%^ = 



AQ%t) < C\a\V%t-) AT%t) < 

as soon as > 4C and 6o < 

When the boundary data changes, then 

AV'ii) < C\Ag'{t)\ AVg'ii) = -\Ag'{i)\ 

AQ'{t) < C\Ag^{t)\V'{t-) AT^{t) < -^\Ag%t)\ 

as soon as Hi > 3C and 6o < ^/H2. 

The above choices are consistent. Indeed, choose first Hi and K, then 
H2 and finally 5o- 

The wave front tracking approximation can be constructed for all times, 
indeed we show that the total number of interaction points is finite: waves 
of families are created only through the Accurate solver and the 

use of the Accurate solver in \^{t,x): t > 0, x > 7(t)} leads to a uniform 
decrease in T*". Therefore only a finite number of waves, which can hit 
the boundary, is present. Since also the jumps in the boundary are finite, 
there are at most a finite number of points in the boundary with outgoing 
waves. This observation, together the argument used in the standard case, 
see [B] , shows that the total number of interactions is finite on all the domain 
{{t,x): t>0,x> 7(t)}. 

As in [6l Paragraph 7.3], the strength of any rarefaction, respectively 
nonphysical, wave is smaller than Ce, respectively Cp. This estimate is 
proved simply substituting Q{t) in [6| formula (7.65)] with the strictly de- 
creasing functional T^(t) defined at ()4.11|) . 

As in the standard case, choosing p sufficiently small we prove that the 
total size of nonphysical waves is bounded by e. To this aim, recall the 
generation order of a wave. Waves created at time t = 0, as well as waves 
originating from jumps in the boundary data, are assigned order 1. When 
two waves interact in the interior x): t > 0, x > 7(t)} of the domain, 
the usual procedure [6, Paragraph 7.3] is followed. When a wave of order k 
hits the boundary, all the reflected wave are assigned the same order k. 

For k >1, define 

jre/j.\ T^ST^ 1 I I has order > k I I , , has order > k I 

V,{t) = ^2. H^^isoffamily<£ ^ ''''^''isof family >£ 
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Qm 




,< , are approaching 
' ' ' one of them has order > k 



As in [6, Paragraph 7.3], for k >1, let 1^ denote the set of those interaction 
times at which the maximal order of the interacting waves is A;. Iq denotes 
the set containing t = and all times at which there is a jump in the 
boundary data. On the other hand, is the set of those interaction times 
at which a wave of order k hits the boundary. A careful examinations of the 
possible interaction yields the following table for k > 3: 







^V§{t) = 
AV^it) < 



t G IoUliU...U/fc_2, 

t G 4_iU4U..., 

t G Ji U J2 U . . . U Jfc_i, 

t G Jfe U Jk+i U . . . . 



Denote the positive, respectively negative, part of a real number by: | x ] 



max{0, x}, respectively \x\ 
mula (7.69)], we get for k>2>: 



-x\ . Therefore, similarly to [6, for- 



0<s<t 



0<s<t 



0<s<t 



Now we need to estimate the last sum: Q|(i) = Y2o<s<t IAQ|(s)l • Ob- 
serve that for k > 3: 



AQlit) + AT%t)-V,%t-) < 

AQlit) + H2 AQl_^{t) ■ V%t-) < 

AQUt) < 

AQUt) + AT%t)-Vi{t-) < 

AQl{t) + AV^it)-V%t-) < 



t G /oU/iU...U4_2, 

t G /fc-i, 

t G /fcU/fc+iU..., 

t G JiU J2U...U Jfe-i, 

t G Jfc U Jfc+i U . . . . 



Hence we can write 



Qm < E 

0<s<t 



-AV,%s)j +H2lAQl_,{s) 



sup V%t) 

0<T<t 



+ E IAT"(^)l_- sup V^ir 



0<s<t 



0<T<t 
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< 35H2 ■ Ql_,{t) . 

By induction we obtain 

QUt) < {?,5H2f-^Ql{t) < {'iSH^f-H. 

Therefore, if 5 is sufficiently small (so that 2>6H2 < 1), there exists A^^ > 
such that the total size of the waves of order greater or equal to A'^^ is smaller 
than e: 

V^{t) < H2 QUt) < H2 {?,5H2f''^ S <e, foi k > . 

Now we observe that the numbers of wave of a given order, is bounded by 
a number that depends on e but not on the threshold p: indeed let the 
maximum total number of waves that can be generated in a solution of a 
Riemann problem inside the domain or at the boundary. Let M be the sum 
of the total number of jumps in the initial data and in the boundary data. 
The wave of first generation are born at i = 0, at the jumps in the boundary 
or when a wave of first generation hit the boundary. Since the waves of first 
generation which can hit the boundary (the ones which belong to the families 
i = 1, . . . ,i) are born only at t = 0, the total number of first generation waves 
is bounded by a constant = M ■ + M ■ ■ not depending on the 
threshold. Suppose now that the number of waves of generation lower or 
equal to k is bounded by a constant not depending on p. The waves of 
order k + \ can be generated only when two waves of lower order interact, or 
when a wave of order k + 1 hit the boundary. Since the waves of order k + 1 
which can hit the boundary can only be generated by interaction of waves 
of lower order, the total number of generation k + 1 waves is bounded by 
C^fc+i = . . Ms + ■ -Ms -Ms which do not depends on p. Fix now k 
such that < e. Hence also the total strength of non physical waves with 
order greater or equal then k is lower than e. Then observe that the total 
number of non physical waves with order less than k is obviously bounded 
by C^. Since the strength of any single non physical wave is bounded by 
Cp, if we choose the threshold p such that Cp ■ < e, we have that the 
total strength of non physical waves is bounded by 2e. Finally we observe 
that if 7 is any ^-non-characteristic curve, then TV (n(-,7(-))) is uniformly 

bounded by a constant time TY{uo)+TV{g). Indeed, this property is proved 
following the techniques in [1'6\ Theorem 14.4.2 and formula (14.5.19)] with 
our strictly decreasing functional T^. 

The following lemma on the regularity of u along non-characteristic 
curves is of use in the sequel. 
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Lemma 4.4 Fix a positive T. Let u he an e- approximate wave front track- 
ing solution to \1.S\) . LetTQ,Ti he l-non- characteristic curves. Then, there 
exists a constant IC > independent from T such that 

[ u{t,To{t)) -u{t,T,it)) < -(TV(no) + TV(5)) ||ri-ro||co([o,T]) 

Proof. Let F be an £-non-characteristic curve. Consider a perturbation 
T] € C'^'''"(M'''; M) with ||7?||co + ||??||loo sufficiently small. By the above con- 
struction of e-solutions, there exist times t^ and states Ua such that 

u {t, T{t)) = ^ X[t^,t^+i[{t) and T{ta) = Kta + Xa (4.12) 



Indeed, here x = A^t + Xa is the equation of a discontinuity line in u crossed 
by r. If {ta,Xct) is a point of interaction in u, then we convene that all states 
attained by u in a neighborhood of [ta^^Xa) appear in the sum in (I4.12p . 
possibly multiplied by the characteristic function of the empty interval. 
If ||?7||qi is sufficiently small, then there exists times t',^ such that 

M (t,r(t) + = ^liaXji^ ^'^d T{t'^) +rj{t'^) = \at'^ + Xa- 

a 

Subtracting term by term, we obtain 

Xa{t'c,-ta) = {r{t'j-r{t^))+7]{t'j 

f « + (1 - di? (4 - 1^) + ?7(4) . 



If' -H 



< 



A, - /; r {M'^ + (1 - ^)ta) dd 
hllc° 



Therefore, 



E 



\u {t, r{t)) - u (t, r{t) + r]{t)) II dt 

— Illf' I ^ ll^llc" \ " 

I I ^ z — J 



< 



TV(n(.,r(.)) < /C^^(TVK) + TV(5)) (4.13) 



proving Lipschitz continuity for rj small. We pass to the general case through 
an interpolation argument. Introduce the map 



T 



u (t, (1 - ^)To{t) + ^Ti{t)) - u (t, T^it)) 



dt. 
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The estimates above prove that the map — > u (•, (1 — '!9)ro(-) + i?ri(-)) is 
continuous in L-*-, hence also ip is continuous and by (j4.13p its upper right 
Dini derivative satisfies 

o.^(,)<cIXM±IXM||r,-r„||,. 

for all ?? G [0.1]- Hence, by the theory of differential inequalities, 

T 

u{t,Tiit)) -u{t,Toit)) dt = V(1)-V'(0) 

TVK) + TV(g) 

completing the proof. □ 

We want now to compare different solutions. Take two e-approximate 
solutions u, V corresponding to the two initial data Uq, Vq and the two 
boundary data g and g. Let a; be a piecewise constant function with the 
following properties: uj{t, •) is an L-'--function with small total variation, 
u{t, x) has finitely many polygonal lines of discontinuity and the slope of 
any discontinuity line is bounded in absolute value by A. The function uj 
does not need to have any relation with the conservation law. 

Define the functions w = v + to and g = (gi, . . . , q„) implicitly by 

w{t, x) = S {q{t, x)) {u{t, x)) 

with S as in (14. 3p . We now consider the functional 

e 



^{u,w){t) = \qi{t,x)\Wi{t,x)dx 



lit) 

"■ r+oo 

+ / \qi{t,x)\Wiit,x)dx 
~e+i J lit) 



(4.14) 



where K is a constant to be defined later and the weights Wi are defined 
setting: 



Wiit, x) = l + x) + K2 (T' (n(t)) + {v{t) 

The functions Ai are defined as follows. Denote by ax^K the size of a jump (in 
u or v) located at x of the family k (k = n+1 for non physical waves). Recall 
that J{u), respectively J{v) denote the sets of all jumps in u, respectively 
in V, for x > 7(t), while J{u), J{v) are the sets of the physical jumps only. 
If the i-th characteristic field is linearly degenerate, we simply define 

A r \ ■ I I I 7/ N , , 7/ N , y < x, i < K < n, or \ 

Mx) = I \ay,.\:yG Jin) UJ{v) and l^^^^^^- 
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On the other hand, if the i-th field is genuinely nonlinear, the definition of 
Ai will contain an additional term, accounting for waves in u and in v of the 
same i-th family: 



Aiix) 



Yl { ky,K| • y ^ -Hu) U J(f) and 



y < X, i < K < n, or 



+ < 



E 
E 



y>x,l<K<i 
y £ J{u), y < X or 



■ y £ J{v), y > X 

y £ J{v), y < X or 
" y £ J{u), y > X 



iiqi{x)<0, (4.15) 
ifqi{x) > 0. 



Recall that non-physical fronts play no role in the definition of Ai. We 
remark that the function u enters the definition of Ai only indirectly by 
influencing the sign of the scalar functions Qi. The constants ki, K2 are the 
same defined in [6j. We also recall that, since 60 is chosen small enough, the 
weights satisfy 1 < Wi{t,x) < 2, hence for a suitable constant C3 > 1, 

^\\w{t)-u{t)\\ <<P{u,w){t) <C3\\w{t)-u{t)\\ (4.16) 

where the L"*" norm is taken in the interval ]7(i), +00 [. 
We state now the following theorem. 



Proposition 4.5 Let the system satisfy the assumptions of Theo- 

rem \2.2[ Then, there exists a constant 6 £ ]0,So[ such that, let u, v, LaJ, w he 
the functions previously defined, satisfying {u{t)), {y{t)), (a;(t)), 
{w{t)^ < 6, for any t > 0, then one has 



<^>{u,w){t2) < <!>{u,w){ti) + Ce{t2-ti) 

"12 



+ C 



b (u (s,7(s))) -b(v (s,7(s))) 



+ TV(a;(s,-)) ds. 



An immediate consequence of the above result that is useful below is 
^{u,w){t2) < ^{u,w){ti) + Ce{t2-ti) 

ft2 



+C j ' {\\g{s) - -g{s)\\ + TV {u{s, •))) ds. 



(4.17) 



Proof of Proposition 14.51 In this proof we use the main results obtained 
in [3l[6]. At each x define the intermediate states Uo{x) = u{x), Ui{x), . . ., 
Unix) = w{x) by setting 



Ui{x) = Si {qi{x)) o S'i_i {qi-i{x)) o ■ ■ ■ o Si (gi(x)) {u{x)) 
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Moreover, call 



Xi{x) = Xi {Ui-i{x),Ui{x)) 



the speed of the z-shock connecting Ui-i{x) with Ui{x). For notational 
convenience, we write gf"^ = qi{y+), (?f~ = qi{y—) and similarly for W^^, 
Af^. If y,y are two consecutive points in J = J{u) U J{v) U J(a;), then 
qf^ = W^'^ = W^~, Af^ = Af~. Therefore, similarly to [6l[l5], outside 
the interaction times we can compute: 



y£j 1=1 

n 



i=l 



Xy) - Wl 



(AZ+ - 7) + E 
1=^+1 



(AJ^ - 7) 



where Xy is the velocity of the discontinuity at the point y. This is because 
the quantities qi vanish outside a compact set. For each jump point y G J 
and every i = 1, . . . ,n, define 



Kq 

9! 



\ii < I 
if i > ^ + 1 

(Ar - %) - ^t" ^Vt - iy) 



so that 



Yl Yl ^V'' 

yeJ i=l 



+KY,wrqr (AZ+-7)+ E 



i=l 



7+ 



i=e+i 



7+ 



(AJ^ - 7) • 



Note that q^^ is a reparametrization of the shock curve equivalent to that 
provided by qf^ and that satisfies the key property, see (Gj Remark 5.4], 

{Si{qi) o Si{-qi)) {u)=u. 

Therefore, the computations in [3, Section 4] and [U Chapter 8] apply. As 
in [21 formula (4.13)] we thus obtain 



ye J i=l 
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Concerning the term on the boundary, (7) imphes that if i < £, then 
A7+ - 7 < -c. Moreover, H^'^+ > 1. Hence, if g"" = b {u{t,-f{t)+) , 
= 6 7(t)+) , Lemma H?2] imphes 



£ n 

< -cKY,\<ir\+cY. 

i=l i=i+l ' 

< -^^E\c\+cj:\<ir 

i=l i=l 



1=1+1 



i^r - 7) 



.7+ 



+ C{U-g^\ + 



,7+ 



provided 

K > C/c (4.18) 
is sufficiently large. Therefore, reinserting the t variable, we obtain 

d_ 
di 



Hu,wm < c(e + TV(u;(t,-)) + |h(t,7W+)|| + ||/(t)-r(t)| 

b(u (s,7(s))) -b(v (s,7(s))) 



< C e + TV(u;(t,-)) + 



Then, standard computations (see [Gj Theorem 8.2]) show that when an 
interaction occurs, the possible increase in Ai{x) is compensated by a de- 
crease in T'=. Therefore, the functional ^ is not increasing at interaction 
times. Hence, integrating the previous inequality, we obtain (I4.17p . □ 



Proposition 4.6 Let system il.^) satisfy the assumptions of Theorem \2.2l 
Then, there exists a process P satisfying [1]) in Theorem \2.2[ 3. in Defini- 
tion \2.1\ and moreover, there exists a positive L such that for all u,v,io, 



\\Pit,to)u-P{t',t'^)v-Uj\\^^ 

< L- <\\u-v-uj\\j^i + \t-t'\ + \to-t'a\ 



+ 



1^'°^' \\g{r)-g{T)\\dT + t-TY{u;)y 



(4.19) 



Proof. Let 5 > be the constant of Proposition 14.51 Define 

Vt = ^ue L^(]R; n) : u{x) = for all x < -f{t) and Tt(u) < 6/2^ . 



20 



Fix Uq € Vq. Approximate the initial and boundary data {uq, Qo) as in (I4.10p . 
Since T^(0) < Tq (u^(0, •)) < To(uo) + e < (5/2 + e < 5, we can construct 
the e-approximate solutions u^{t,x). As in j6l Section 8.3] we observe that 
for < e' < e, the e'-approximate solution is also an e-approximate solution. 
Therefore, we can apply (14.170 with in place of v with lo = and g = g. 
Hence, because of (I4.16|) . we obtain 



u'{t) - (t) 



< L ■ 



+ e-t. 



For any t > 0, u^{t) is a Cauchy sequence which converges to a function 
u{t) E Li(M;M") that vanishes for x < 7(t). The potential T%t) defined 
on e-approximate solutions is non increasing and differs from due 
to nonphysical waves and to the different boundary conditions g^ and g. 
Therefore, ()4.10p and the lower semicontinuity of the total variation and of 

implies that u{t) £ T>t for any t > 0. We set 0)mo = u{t). It is 
obvious that our procedure can start at any time to > 0, so we can define 
P{t,to)u G Vt+to for any uGVt^. 

We want to show now that the map just defined satisfies all the prop- 
erties of Theorem 12.21 The Lipschitz continuity t — s- P{t, to)u is satisfied 
by construction. If we now consider a different initial and boundary data, 
say {v,g) and the same boundary curve 7, in general we have a different 
map P. Taking the limit in (j4.17|) and using (|4.16p for the corresponding 
e-approximations, we get that for any L"*" function uj dependent only on x 
and with small total variation 



< L 



P{t,to)u- P{t,to)v - UJ 
\\u — V — Uj\\i^i + 



l'''^'\\g{T)-g{T)\\dT + t-TY{uj)^ 



(4.20) 



bounding the dependence from the error term uj and proving the Lipschitz 
continuity in g and u. 

Point 3. in Definition 12.11 is obtained by standard methods, see (Gj Sec- 
tion 7.4]. 

Concerning the process property, take u G "Dq and consider its e-approx- 
imation u^. Let e € ]0,e[ and call be the e-approximate solution with 
initial datum u^{t) at time t. Then, if s > 0, is also an e-approximate 
solution in [t,t + s]. Therefore, applying ()4.16p and (I4.17P in the interval 
[t, t -|- s], we obtain 



P{t + s, 0)u - P{s, t) o P{t, 0)u|| 
: + s)-P(s,tK(t)||Li 

u%t + s) - u^{t + s) 



lim lim 
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< lim lim C 
= 0. 



t+s 



+ e s 



We can repeat the same argument for any initial data to > 0. 

Concerning the dependence on the initial time to, take < to < ^'^d 
u G Vt^, u' G Vti^. If < t < tg — to, then obviously 



\\P{t,to)u - P{tX)u'\\^^ < C [\to~t'^\ + 
If t > tg — to, the process property implies 

\\P{t,to)u- P{t,t'o)u\\-^^, 
= ||P(t + to - t;, t'o) o P(t; - to, to)u - P{t, t'o)n'||Li 

< C ||P(t;- to, to)^^ - ^'IIli +C|to-t^| 

< C\\u-u'\\^^+C\to-t'^\, 
completing the proof of ()4.19p . 



□ 



The following proposition extends to the present case the key properties 
of the Glimm functional (14.91). 



Proposition 4.7 Let system satisfy the assumptions of Theorem[ 
Then, for any u G Vq, the map t ^ Tt (P(t, 0)u) is nan increasing for 
t > 0. 

Proof. Above, we showed that the map t T^(t) decreases along e- 
approximate solutions. The monotonicity of t ^ [P{t, 0)m) follow pass- 
ing to the limit e — > 0, thanks to the lower semicontinuity proved in [11|9], 
to (|4.1U|) and to the lower semicontinuity of the total variation. □ 



In order to complete the proof of Theorem [221 we prove propositions 12.3 
and 12.41 together with an auxiliary lemma. 

Proof of Proposition [273l Let be an e-approximate wave front tracking 



solution converging to u. Since the convergence is also in Lj-^^(I 



'7) ■ 



apply 



Lemma 14.41 and Lebesgue Dominated convergence Theorem to obtain: 
rT 



dt 



u{t,To{t)) -u{t,Ti{t)) 
lim^ / / u(t,ro(t) + x) -u(t,ri(t) + 

{t,ro{t) + x) - (t,ri(t) + 



dt dx 



S rT 



lim lim -- , , 

(5-^0 e-»0 Jq Jq 



X 



dt dx 
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< ;c.TVK) + TVfa) 

C 

completing the proof. 



Ollco([0,T]) 



□ 



Lemma 4.8 Letu^ be an e- approximate wave front tracking solution to U.^) 
converging to u. Let T be an i -non- characteristic curve. Then, 



u 



u 



,r(-)) 



-'locV 



■n). 



Proof. By the convergence of to u in Lj-^p(D^; there exists a sequence 
Ejy converging to such that for a.e. x 



u 



{■,T{-)+x)^u{-,,V{-)+x) 



m 



-'loci 



Then, for any T > and for any x for which the convergence above holds, 

cT 

' " dt 



< 



u'- {t,T{t)) -u{t,r{t)) 

u^" {t, T{t)) - u^" {t, T{t) + x) dt 

+ [ u'- {t, r{t) + x) -u (t, r{t) + x) 

Jo 

+ I u{t,r{t) + x) -u{t,T{t)) 
Jo 



dt 



< 2clXM±IM|,| + 



T 



U 



C JO 

where we used Lemma 14.41 and Proposition 12.31 Hence 



t,T{t)+x) -u{t,T{t)+: 



dt 



lim sup 

u^+oo Jo 



u'" {t,T{t)) - u {t,T{t)) 



dt < C\x\ 



and the final estimate follows by the arbitrariness of x, independently from 
the sequence £iy, thanks to the uniqueness of the limit u. □ 



Proof of Proposition 12.41 Let u*", respectively n^, be an e-approximate 
wave front tracking solutions of (|1.2|) . respectively (|1.3p . Apply Proposi- 
tion 14.51 and use the equivalence (I4.16|) to obtain 



7(i) 



\u^{t, x) — u^{t, x)|| dx 



< L 



+ 



7(0) 



|n'^(0, x) — u^(0, dx 



ds + Cet 



and the limit e ^ completes the proof. 



□ 
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Proof of Theorem 12. 2L To conclude the proof of Theorem 12.2^ observe 
first that for the e-approximate solutions, we have 



< £ 



therefore as e — > Lemma 14.81 implies [^). 

Finally denote by let P^T'S) and v'l''^^ the process and the domains 
corresponding to the boundary curve and data (7,3). Fix two boundary 
curve and data (7,5), (7,5), two initial data Uq G 'D^'^\ Uq € V^q''^^ and 
define 



min{7(t),7(t)}, Ti{t) 
jo for x< Ti{0) 

] Uo{x) forx>ri(0) 



ro(t) 

Uo{x) 

m 

By Proposition 12.41 we have for x > Ti{t): 



max{7(t),7(t)} 

Jo for X < ri(o) 
I Uo{x) for X > ri(0) 

p(^'SHt,o)uo\ (ri(t))). 



'p^'''S\t,0)uo 


(x) = 


' P^^^~^\t,Q)Uo 


(x 


P^^~^\t,Q)uo 


(x) = 




(x 



Applying the result for the unchanged boundary curve, we get: 

P^'<^3\t,Q)Uo-P^^~^\t,Q)Uo 



r()(t) 



dx 



+ 



ri(t) 



dx 



P^^^''3\t,0)Uo - P'^^'~^Ht,0)^o 

< C\ri{t) - To{t)\ + C f \\~g{t) - f(t)|| dt + C\\uo - ^o||li 

JO 

< CIlFi - Follco + C||no - ■UoIIli 
\m - git) + g{t) - -g{t) + m - f dt 



+C 



< (^117-71100 + G|I^^o-^^o|Ili + C / \\9{t)-git)\\dt 



+C 



~g{t)-9(.t) dt + C / g{t)-g{t) dt 



Finally Proposition 12.31 and the Lipschitz continuity of b imply 



Igit) - 9it)\\ dt 
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p(^'5)(t,o) J (ri(t)) ]-b( [p(^'^)(t,o)nol (^W) 



dt 



< c||ri-7|lco< 



I CO 



completing the proof of since the computations for g and g are identical. 

We prove now the tangency condition^j. Fix to > 0, u G T>t^ and let F be 
defined by ()4.6p and denote by P the process defined above with g replaced 

o + i)) ) with u as in ()4.5p . By Proposition [2 



by g{to + t)=b 

F{t,to)u = P{t,ta)u. Using 0j, we have 



t 
L 

< — 

~ t 



\\P{t,to)u- F{t,to] 



u\ 



to+t 



to 

rto+t rto+t 

-J \\9(.s)-g{to+)\\ds + j / 

J to ^ to 



\9{s) - g{s)\\ ds 



\9{to+) - g{s)\\ ds. 



The latter term vanishes as t — s- by the definition oi g{to+)- Consider now 
the former term. Fix a positive and sufficiently small 5 so that the curve 
o) is non-characteristic. Let ^ G [0;1]- 



7 / Wais) -g{to+)\\ ds 

J to 

6((5,_<» (7(s))) -g{to+) 



< 



ds 



to+t 



to 



to 



b [{Ss-tou) (7(s))) - b [{Ss-tou) ((1 - ^Ms) + ms))) 

b [{Ss^tou) ((1 - Chis) + - 9{to+) 



ds 



By Proposition 12.31 the first term is bounded by 



C 



7- ((l-OT + eV') 



C 

C0([t„,t.+t])-7"^"^"c''([*'"*°+*l) 



< C6. 



Concerning the latter term, integrate on ^ over [0, 1] and obtain, with the 
change of variable x = (1 — C)7('S) + and with u'^ as in Lemma |4.H 



to+t 



to 



b [(Ss^toU) ((1 - Olis) + ms))) - 9{to+) 



to+t 



i Jto V'(S) - lis) J-y(s) 



6((S,_t„n)(x))-6K: 



ds 
dx ds 



C 

< — 
- tb 



to + t 



to 



S-ta 



O J-f{s) 



I (S's_t„n)(x) — ti'^ll dx ds 



(4.21) 
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Following [HI Section 9.3], let C/^ be the Lax solution to the Riemann problem 

dtu + d^f{u) = o 

u{0, x) = 



if X < 
u{-i{to)+) ifx > 0. 



By the basic properties of the solutions to Riemann problem and the defi- 
nition of V', for all s G [toit] and x € [■j{s),'ip{s)], U'^{s,x) = . Then, 



C 



to+t 



1 



J to S - to J^(s) 



{S,-toU)ix)-U\s,x) 



dx ds 



l-to+t 

By [6l formula (9.16)], 



< ^ 
- t6 



-i(to) + {s-to)\ 



7(to)-(s-to)A 



(S,_t»(x)-C/«(s,x) 



dx ds 



lim 



1 



l{to) + ('^-to)\ 



{Ss-toU){x)-0{s,x) 



S^to S-to J^(to)-(s-to)X 

so that lim^^n ()4.2ip = 0. Collecting the various terms, 

rto+t 



dx = 



lim sup 



iup- / \\g{s) - g{to+)\\ds < C6 

t Jto 



and by the arbitrariness of 6, the tangency condition follows. 

The characterization of P through [7p, \^ with u; = and implies its 
uniqueness through standard computations, see for instance (Sj Section 6, 
Corollary 1]. □ 



5 The Source Term 

This section is devoted to the source term, similarly to O [71 El [HI [TH] but 
following the general metric space technique in [10], applied to L-*- equipped 
with the L-'--distance d. The key point is to show that the map 

F{t, to)u = P{t, to)u + t G {Pit, to)u) X[^(,^^,)^^^[ (5.1) 

is a local flow in the sense of [101 Definition 2.1] on suitable domains and 
satisfies the assumptions of [lUl Theorem 2.6]. 

Following [9, Section 3], we modify the functional $ in ()4.14p and define 
$t on all piecewise constant functions, not necessarily e-approximate solu- 
tions. Therefore, the definition of does not consider nonphysical waves 
and $0 = $ at time t = 0. Consider two piecewise constant functions 
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u,v {]^{t), +oo[, M") with finitely many jumps and assume that TV(u) 
is sufficiently small. 

Define q = {qi, . . . , qn) implicitly by 

v{x) = S {q{x)) (n(x)) 

with S as in (14.31). We now consider the functional 



p+oo r+oo 

^t{u,v)=KS^i \qi{x)\W i{x) dx + "S^ / \qi{x)\W i{x) dx 

where K is defined in the proof of Proposition 14.51 and the weights Wi are 
defined setting: 

Wiix) = 1 + KiAiix) + K2 (Tt(n) + Ttiv)) . 

The functions Aj are defined as follows. Let ax,n be the strength of the K-th 
wave in the solution of the Riemann problem for ()4.7p in u or w located at x 
of the family k. Differently from the notation in SectionfH J{u), respectively 
J{v) denote the sets of all jumps in u, respectively in v, for x > j{t). Indeed, 

we let X = 7(t) in J as soon as b (^u (7(i)+)) 9{t) ^'^'^ the waves cr^^^i-j^j^ 
are defined as in (j4.8|) . 

If the i-th characteristic field is linearly degenerate, we simply define 



Ai{x) = ^ S \^y,K\- y e J{u)U J{v) and 



y < X, i < K < n, or 
y>x,l<K<i 



On the other hand, if the i-th field is genuinely nonlinear, the definition of 
Ai will contain an additional term, accounting for waves in u and in v of 
the same i-th family: 

Ell I y £ Jiu), y < X or \ -r / \ 

I ' ^' ' ye J[v), y > X I 



+ < 



Ell I 2/ £ Jiv), y < X or / % ^ „ 

I I ^'1 y G J{u), y > X I '^n ; - 



The constants ni, H2 are the same defined in [6l Chapter 8]. We also recall 
that, since 6o is chosen small enough, the weights satisfy 1 < Wi{x) < 2, 
hence for a suitable constant C3 > 1 we have 

7^ \\v - njlLi < ^t{u,v) < C3 \\v - -uIIli, 
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where the L"*" norm is taken in the interval ]7(t), +00 [. 
For a fixed positive M, define 



M 



u{x) = for all x < 7(t) 
uGL^(M;f7): rt{u)<6-C(T-t) 
\\u\Ui < Me^^ + Ct 



with tt defined in (|4.9p . C, 5 and T to be specified below. 

Lemma 5.1 For all to G [0,r], t > sufficiently small and u,u (z T>t^, 

rt^+t[F{t,to)u^ < TtM + ct 

^to+t(F{t,to)u,F{t,to)u) < {l + Ct)^tSu,u). 



(5.2) 

The proof is as that of [7, Lemma 3.6 and Corollary 3.7], see also [HI 
Lemma 2.3]. 

Corollary 5.2 For t small, F in lO]) satisfies F{t,to)'D^^ C P^jf+i- 

Proof. The bound on is a direct consequence of Lemma f5.ll Concerning 
the estimate on the L"*" norm, for u € i>f^ , compute: 

F{t,to)u 

L 

< \\P{t,to)u-u\\^^ + \\u\\J^^+t\|^G{P{t,to)u) -G{0) +i||G(0)||Li 

< Ct + \\u\\j^i + Ct {\\u\\-^i + Ct) + Ct 

< {1 + Ct)\\u\\^^ + Ct 

< (l + Ct) (Me^*° + Cto) + Ct 



< Me^^*''+*'> + C{to + t) 
hence F{t, to)u is in 



□ 



In what follows, relying on [10^ Condition (D)], we consider F as defined 
on the domains I?^^ and not on a single domain, as in |101 Definition 2.1]. 

Proposition 5.3 The map F defined in \5.1\) is L-*- Lipschitz continuous, 

satisfies F{0,to)u = u for any {to,u) G |(r, ru): r € [0, T], w € and 

there exist positive C, independent from M, such that for to,t'^ £ [^jF], 
te[0,T- to], t' G [0, T - t'^], u G V^^, v! G Vf 



F{t', t'o)u' - F{t, to)u ^ < L \\u - nll^i + (l + lA^-^) \t! - t| + jt^ - to\ 
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Proof. Compute: 



F{t',Ou' - F{t,to)u 



< \\Pit', Qu' - P{t, to)u\\^^ + \t' - t\ 



G(P{t',t')u')Xr , , , 



+t 



G {Pit'X)u') - G {Pit,Qu) Xj^,^,,,_,^j 



< \\P{t',t'^)u' -P{t,to)u\\^, 
+ \t' - t 

+t 



+t 



G{P{t',Cu-G{0)) ^^ + ||G(0)||l. 
G{P{t',Ou') -G{P{t,to)u " 



),+oo[ -''^[7{to+t),+oo| 



< (1 + Ct) \\P{t',t',)u' - P{t,to)u\\^^ + \t' - t\ ||G(0)||li 
+C\t' -t\ \\P{t',t'^)u\\^^ +Ct\j{t'^ + t') -j{to + t)\ 



< C y\t' — t\ + \t'g — to\ + — 

+C\t' -t\ {\\P(t',Qu' -u'\\^^ + \\u'\\^^ 



< C\^[i + j \t' -t\ + \t'^ - to\ + \\u' - 

completing the proof. 

Recall \10\ Definition 2.3]: an Euler e-polygonal is 



n 



k—l 

F%t, to) u = F{t- ke, to + ke)o Q F{e, to + he) u (5.3) 

h=0 

for k = [t/e]. Above and in what follows, we denote the recursive composi- 
tion Or=i fi = fi ° f2 ° ■ ■ ■ ° fn- Here, [ • ] stands for the integer part, i.e. for 
s € M, [s] = max{fc £Z: k < s}. 

The hypotheses to apply [lOl Theorem 2.6] are satisfied. 



Proposition 5.4 The local flow F in 15. 1\) is such that there exist 



M 



1. a positive constant C such that for all to G [0,T] and all u G 
d {F{kT, to + T)o F{t, to)u, F{{k + l)r, to) <CkT'^ 
whenever k gN, {k + 1)t, r G [0, T — to]; 
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2. a positive constant L such that 

d (^F%t, to)u, F^{t, to)u?j < L ■ d{u, w) 

whenever e e]0,5], u,w e Vf^ , t > and to, to + t e [0, T]. 

Note thatd states that [10^ 1. in Theorem 2.6] is satisfied with uj{t) = Ct. 

Proof. To proved, the key property is 0j in Theorem 12.21 see also [TT| 
Proposition 4.9]. 

F{kT, to + t) o F(r, to)u -F{{k + l)r, to) u 

= P{kT, to + r) {p{t, to)u + tG {P{t, to)u) X[^(^t.+r),+oo\ 

+kr G (^P{kT, to + r) (^P{t, to)u + tG {P{t, to)u) X[^(t„+,),+oo 

■^[7(to + (fc + l)T), + Oo[ 

-P{{k + l)T,to)u 

-{k + l)r G (P ((A; + 1)t, to) n) X[,(^t.+^k+i)r),+oo[ 
= P{kT, to + r) ( P(r, to)u + tG (P(r, to)u 



7(to+T), + Oo[ 

Pikr, to + r)o P{r, to)u - r G (p {{k + 1)t, to) u) 



+kT 



G P{kT, to + t) ( P(t, to)u + tG {P{t, t. 



-G{P{kT,to + T)oP{T, to 



X[7(io+{fc+l)T),+oo[ • 



Using in Theorem 12.21 in the first two hnes with t = t' = kr, to = t'^ 
for to + T, V = P{T,to)u, to = tGx and in the latter two lines (G), in 
Theorem 12.21 with uj = 0. We thus get 



d (F{kT, to + T)o F{t, to)u, F {{k + l)r, to 
F{kT, to + r)o F{t, to)u -F{{k + l)r, 



I u 



< Gt 



G{P{T,to] 



•o+t),+oo[ 



G(P{{k + l)T,to)n]xj,u 



7{*o + (fe+l)r),+oo 



+CkT^ 



G{P{T,to, 



u) 



LI 
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< Ct {kT ||G||loo IItIIl- +CkT)+C (||G(0)||li + 1 + m) kr"^ 

< C{l + M)kT^. 

The bound [B is a direct consequence of the equivalence (14.160 and (15. 2p in 
Lemma |5.H see also |1H Proposition 4.9] and fP; formula (3.1)]. □ 

Proof of[^, \^ and fl^) in Theorem 13.21 By [ini Theorem 2.5], for any 
M, the local flow F generates a Lipschitz process P on 15*^. By the char- 
acterization of P as limit of Euler polygonals, it follows that P is uniquely 
defined on all 

a/^o rt{u)<6-C{T-t) f- 

Hence, P satisfies [Hj in Theorem 13.21 and Fz)) holds. 

To prove [i^j, note that j F{t.to)u — F{t, to)u — as t ^ 0, for all 

u € Vt^ and apply [10, cj in Theorem 2.5]. □ 
For any iV G N, define the operator Hat: L1(R;M") PC(R;R") by 

-1+^^ j-(k+l)/N 

Lemma 5.5 IlAr is a linear operator with norm 1. Moreover, TV(nArn) < 
2TV(m) and for all u G L1(M; M'^) n BV(M; M"), Unu ^ u in V-{R;W'). 

For the proof, see [3 Lemma 3.4]. 

Proof of Proposition 13.31 Set for simplicity to = 0. Let e, e > and 
A'^ G N be fixed. Consider an e-approximate wave front tracking solution 
ye,e,N _ ^e,e,7V^^^^^ (jl.2|) on the time interval [0,e[. Define it at time 
t = e setting 



Extend n^'^'^ recursively on [0,T]. Note that 

lim lim u''^'^(t)= F^(t,Q)uo 

where F^ is defined in (|5.3p and Uq is the initial datum in (jl.ip . Note that 
this is the usual operator splitting algorithm. 
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Given any curve -^-non-characteristic curve T with support in D^, define 



n+l 



\ a;>r(t) i=l 7{t)< 

+TV (n^'^'^(-,r(.));[0,t] 



(5.4) 



for suitable positive K, K. 

Computations similar to those above allow to prove that < 
) for all t ^ eN. Indeed, when a wave crosses F, the increase in 
TV(u^''^'") is compensated by the decrease in the first term on the right 
hand side of ()5.4p . 

At times t G eN, < H^(t-) + Ce. Therefore, for all t G [0,r], 

< H'^(O) + Ct. By (|5.4|) . we get that there exists a C dependent only 
on Uo and i such that 

TV U'''^{;T{.));[0,t])<C. 



Consider now two ^—non-characteristic curves ri,r2 with support in B^. 
The same steps in the proof of Lemma 14.41 lead to 



T 



U 



e,e,N 



u 



e,e,N 



c 

dt< — ■ iir 



Til 



CO 



Let now e ^ and N +oo, with the same technique of Proposition 
we obtain 



T 



F'{t,0)u (ri(t)+) - F'it,0)u {T2{t)+) 



dt < ^IlLa -Lillco (5.5) 



with a constant C that now depends also on T and on L2 in (G). Let now 
also e — > and, by the LJ-^^(]R^; $7) convergence of the Euler polygonals, 
obtain as in Proposition 12.31 that 



T 



P(t,0)nj (ri(t)+) - [P{t,0)u) {T2{t)+) 
completing the proof. 



(J 

dt < — IlLa -r 
c 



ilico (5.6) 
□ 



Proof of Imp and[^ in Theorem 13.21 The Lipschitz continuity upon the 
initial data is a consequence of [101 , b) in Theorem 2.5], thanks to Proposi- 
tion 15.41 The dependence of the Lipschitz constant for the variable t on the 
L-*- norm of the initial data is shown in Proposition 15.31 

The Lipschitz conditions (|5.5p and (j5.6p allow to prove the L-*- ([0, T]; J7) 
convergence of the traces as in Lemma 14.81 

[f%;0)u) (r(.)+) - (P(-,0)n) (r(.)+) in L\[0,T];Q). 
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The map t F^(t, 0)u satisfies for a.e. t the boundary condition, hence the 
same does the solution t P{t, 0)u, proving 2. in Definition 13.11 Condi- 
tion 3. in the same definition is proved using the tangency condition [3|) , as 
in [3 Corollary 3.14]. 

We are left to prove the Lipschitz dependence from the boundary and 
the boundary data. To this aim, introduce two boundaries 7 and 7, with 
7 < 7 and boundary data g, g. Let Vt, 'Dt, P^''^{t-,to) and P^'^{t,to) the 
corresponding domains and processes. We need to prove that for any u E 
Vq n Vq (and therefore u{x) = for x < 7(0)): 



P^''^(t,0)n-P5'^(t,0)M 

Note first that 

P3'^(t,0)'u-P''^(i,0)'u 



< C 



7-7llco([o,t]) + / \\9{r) - g{T)\\dT 



co([o,t]) 



+ 



P^'^(t,0)M-P^'^(t,0)u 



where It = [7(t),+oo[. Hence, we consider below only the latter term in 
the right hand side above. Introduce the linear projector vr^f = vXr 

denote w{t) = P^''^{t,0)u. Then, applying [6l Theorem 2.9] to the process 
P^'^ and to the Lipschitz curve r — > tTt-w{t), using the tangency condition, 
we compute 



P<^'^{t,0)u- pS'^it,0)u 



< L I liminf 



< L I liminf 

£^0 



< L I liminf 



7rr+eW{T + e) - P3^^{e, t) {TlrW{T)) 



P9^"l{e,T)w{T) - P9^^{e,T) {TTrW{T)) 



dr 



P3'^e,T) {w{t)) - pS'^e,T) {7:Mr)) 



dT 



dT 



+ L 



G {w{t)) - G {-K^wiT)) 



e 

dT. 



For the term deriving from the source, we use the L"*" Lipschitz continuity 
of G to estimate: 



G {w{t)) - G {'k^w{t)) 



dT < C 



w{t) — (7rT-ii;(r)) 



dT 



< Cr||7-7||f.o 



c°{[o,t])- 
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Concerning the other term, denote by F^°'"'{t, to)u the tangent vector defined 
in (j4.6p . Here, we expHcitly denote the dependence of the tangent vector 
on the curve 7 and on the pointwise boundary data Qo = g{to)- By [3]) in 
Theorem 12.21 the curve rj P^'^{'q,T) (7rT-t(;(r)) is first order tangent to 
•q — > F^^'^^'^ (r],T) (7rT-w(r)), while rj P^'^ {ri,T)w{T) is first order tangent 
to ry — > F^^'^')'"'{r],T) (u;(r)). Because of the finite propagation speed, the 
two tangent vectors coincide in the interval [7(t) + r/A, +00 [. Therefore, 



Jo 



lim inf 



F9W'T(e,T)u'(r) - F9W'^(e,T) {-Krw{T)) 



dT 



lim inf 



t n{r)+eX 

lim inf — / 

e J^{T+e) 



L^(/r + .) 



dT 



F5M'T(e,r)ii;(T)) (x) 



dxdr . 



Referring to Lemma 14.11 introduce the quantities 



w'^{x) 
w'^ix) 



w (t,7(t)) , 

Ipni^n) O . . . O '0f+l(<T£+i)(u;'^'^). 

{w{t, x) for X > 7(t), 

t(;(r, 7(r)) for 2; < 7(t), 

{w{t,x) for X > 7(r), 
^<J,T £qj. ^ ^ ^^^^ 

By formulae ()4.5p - ()4.6p and since the boundary condition is satisfied for 
almost all r, that is b (^w (t, 7(r))^ = ^(t), one has for x > 7(t + e) 

(FfW'^(e,r)u;(T)) (x) = (SsW^ix) 

(F5W'^(e,r)(7r,t«(r)))(x) = {Sew){x) 

where 5 is the purely convective Standard Riemann Semigroup without 
boundary generated by / [S| Definition 9.1]. 

Denote by C/| and C/| the solutions to the two Riemann problems: 

Ut + f{u)x = 

w'^{^{t) — ) for X < 
{li'^)) for X > 



n(0, x) 



Ut + f{u)x = 

|tS^(7(r)-; 
(7(r)) 



ti(0, x) 



for X < 
for X > 



34 



Formula [HI (9.16)] implies that 



/ lim inf — / 



{S.W^ (X) - {SeWn (X) 



dxdr 



/■* 1 

< / lim inf — 

Jo e J7(T+e)-7(T) 



C/|(e,x)-^|(e,x) 



dxdr 



By Remark 13. 41 for almost all r such that 7(t) < 7(t) one has w'^ {'j{t)—) = 
w'^{^{t)—) = wl, therefore C/|(e,x) = Wr- While for almost all r such that 
7(r) = 7(r), the boundary condition implies w'^ {^{t)—) = w'^ (7('7")) = 'u^r 
therefore we have again C/|(e,x) = wl- We compute, for almost all r 



C/|(e,x)-C7|(e,x) 



< C7 



C 



E^{w;,g{r))-E^{w;,b{w;)) 



< c 



g{T) - b (u;;) 



Finally we compute, using Proposition 13.31 

[/|(e,x)-C7|(e,x) 



lim inf — 



£-♦0 e 



7(r+£)-7(r) 



dxdr 



< C 



g{T) - b {wj.) 

< C 1^ \\g{T) - g{T)\\ dr + cj^ b (w {t,j{t))) - b (w {t,j{t))) 



dr 



<C ||fi((T) -5((t)|| (iT+ ||7- 711^0. 
Jo 

The general case of two non ordered curves follows immediately by the tri- 
angle inequality. □ 
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